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ABSTRACT 

Dipole representations of the earth* s magnetic 
field have been found to have insufficient accuracy 
for the study of magnetically trapped particles. A 
coordinate system consisting of the magnitude of the 
magnetic field B, and the integral invariant I has 
been found to adequately organize measurements made 
at different geographic locations. It is shown in 
the present paper, that a parameter L = f(B,l) can 
be defined which retains most of the desirable 
properties of I and which has the additional property 
of organizing measurements along lines of force. Since 
the parameter L is the analog of a physical distance 
in a dipole field (the equatorial radius of a magnetic 
shell), it is usually found to present fewer conceptual 
difficulties than the integral invariant I. 
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INTRODUCTION 

Considerable difficulties were encountered in 
the early attempts to map the measured intensities 
of trapped particles. The difficulties were chiefly 
due to the high spatial gradients in the intensities. 
At low altitudes in the inner zone for example, a 
change in the distance from the center of the earth 
of 3 percent results in an intensity change of over 
a factor of 10. The shape of the earth* s field 
deviates more than 5 percent from any possible dipole 
representation, therefore attempts to use a dipole 

model for a coordinate system can be expected to 

> 

produce rather chaotic results. 

The present paper describes a coordinate 
system which not only takes account of the nondipole 
character of the field but can also be used to organ- 
ize measurements along lines of force as is necessary 
for most theoretical studies. 
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ADIABATIC INVARIANTS 


The use of adiabatic invariants in describing 
the motion of a particle in a magnetic field has been 
discussed by Northrup and Teller /19527* In the 
present paper, the integral (or longitudinal) in- 
variant at a point in space A is defined as 
A* i /o 

I - <p (l B 0 /B) '^ds where ds is the differential 
A JL 


path length along the line of force connecting the 
point A with. its conjugate point A*, B ^ is the magni- 
tude of the magnetic field along the line of force, 

B is the magnitude of the magnetic field at the point 
A, and where the integral is to be taken along the 
line of force between the conjugate points. Defined 
in this manner, I can be considered to be a scalar 
field which has a definite value at each point in 
space and does not require reference to the motion 
of trapped particles. If there is no electric field. 


this definition is equivalent to I = P 1 P ds 

A « 


A* 


where P is the total momentum of a particle which 
mirrors at the point in space A and where P ^ is the 
component of momentum along the line of force. If 
we assume that all three adiabatic invariants are 
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conserved for a given particle, then the values of 
I are the same at every point in space at which a 
given particle mirrors. If the energy and mirror 
point distributions of a set of trapped particles in 
time equilibrium do not change importantly during the 
time the slowest particles drift once around the 
earth (for example, the high energy protons in the 
inner zone), then the directional intensity perpen- 
dicular to the line of force will be the same at 
all points in space which have the same values of B 


and I. 



6 


MAGNETIC SHELLS 

The points in space which have the same value . 
of B and I form a ring in each hemisphere. A particle 
mirroring at this B and I will remain upon the surface 
(or "shell") described by the lines of force which 
connect these rings. 

Now in general, two particles which initially 
mirror at different values of B along a particular 
line of force will not drift in longitude to the same 
lines of force. This means that the shells described 
by the trajectories of the two particles do not 
coincide. It will be shown however, that in the 
earth's magnetic field, this effect is relatively 
small. One important consequence of this fact is 
that the omnidirectional intensity as well as the 
directional intensity. is constant along the loci of 


constant B and I. 
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DEFINITION OF THE MAGNETIC SHELL PARAMETER L 

The fact that all particles which drift through 
a given line of force will remain on approximately the 
same shell throughout their motion leads immediately 
to the desirability of finding a method of labelling 
all points in space with a number which is unique for 
each shell. 

Each locus of constant B and I is a line which 
lies on a fixed shell. Any parameter which is a 
function of only B and I will therefore have the correct 
longitude dependence. The problem of defining a mag- 
netic shell parameter is thus reduced to finding a 
function f(B,l) which is constant along lines of 
force. Such a function could be found empirically by 
examining the functions 1^ = g^(B) along lines of 
force calculated with a representation of the earth* s 
field. Ah empirical fit is not necessary however, 
because the average of the functions g^(B) around a 
shell is well represented by the function obtained 
for a dipole field. 

The integral invariant in a dipole field at 
a magnetic latitude A on a line of force which has an 
equatorial radial distance ol R Q is given by 
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I = 


or 



1-Y 2 




1/2 


, P V 2 

)(l+3Y a 2 ) dY 


a 

(i) 


i = R 0 h i( x ) ( 2 ) 

where Y = sin A. Now in a • dipole field 

1/2 

B = ■ -3 M 6 (1 + 5 Bln 2 X) (5) 

R 0 P ooa D A 

where M is the dipole magnetic moment. By (j3) we 
see that 

R 0 Vm = h 2 (A) (4) 

Equations (2) and (4) give 

A/M = R 0 ^B/M h 1 ^(A) = R 0 5 B/M h 5 (R 0 5 B/M) 

or 

r 0 5 b/m = f( A/m) (5) 

The magnetic shell parameter L is now defined 
for a point in the earth's magnetic field by the 
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equation 

L 5 B/M = F(I 5 B/M) (6) 

where I and B are to be calculated for the point 
with a representation of the earth's field, M is 

-ft 

the dipole moment of the earth and F is the function 
in equation (5) which is calculated with a dipole . 
field. 

A set of values for the function F is given 
in Table I. For accurate computation of L, the 
following method can be used: 

Let 

ln(^ - 1) = a n X n (7) 

where X = ln(l‘ 5 B/M). Sets of the coefficients a n 
for different ranges of X are given in Table II. 

This method introduces an error in L of less than 
0.3 percent for -co<x<co and less than 0.03 percent 
for X < 10. 

* The value M = 8.06 x 10 2 ^ gauss cnr 5 = 0.311653 gauss 

R ^(R = 6371.2 km) has been found to be satisfactory. 

6 0 
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SIGNIFICANCE OF L 

In general, L should be regarded as a parameter 
which retains most of the useful properties of I 
exactly and which is also approximately constant along 
lines of force. To gain some feeling for the meaning 
of L, it is useful to remember that in the case of a 
pure dipole field, a magnetic shell is labelled by 
an L equal to the shell's equatorial radius. 

For most purposes, B and L should be used as 
the principal spatial coordinates. For studies 
involving intuition, it is sometimes desirable to 
use a coordinate system which has a closer resemblance 
to the actual physical geometry. One system of this 
kind can be obtained by transforming B and L to polar 
coordinates using the dipole relations 

B . J*. U - fj /S , R = l cos 2 (8) 

This relationship is illustrated in Figure 1. 

A great deal of caution must be exercised 
when using the coordinates R and 7 \ obtained in this 
manner, because while the irregular characteristics 
of the magnetic field will in effect be removed, 
geographic coordinates will transform in an irregular 
longitude dependent manner. 
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VARIATION OP L ALONG LINES OP FORCE 

Tables of I versus B along 1400 different lines 
of force were kindly supplied to the author by 
Messrs. D. C. Jensen, R. W. Murray, and J. A. Welch, Jr. 
These tables were computed with the 512 term spherical 
harmonic expansion (n = 24, m . = 17) of the earth»s 
field obtained by D. C. Jensen and W. A. Whitaker, 
/19607 and from the 1955 surface data. The L values 
obtained from these sets of B and I were found to 
vary less than one percent along most of the lines 
of force. Figure 2 shows the variation in L along 
six lines of force and includes several of the worst 
cases. The variations along different lines of 
force were found to be systematic in both latitude and 
longitude. This is illustrated in Figure 5 where the 
deviations in L, along, lines forming the magnetic 
shells labelled by L = 1.2, 1.5, 2.2 and 5*0 earth 
radii, are plotted versus the geographic longitude 
at the northern ends of the lines. A positive 
deviation means that L increased toward higher values 
of B. 

This kind of check upon the variation of L 
along lines of force tends to be independent of the 
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precision of the magnetic field analysis. An 
analysis of the earth's field based upon improved 
knowledge of the surface field would probably have 
nondipole terms of similar magnitude and therefore 
would be expected to give similar variations in L 
along lines of force. It should be noted however 
that in the present analysis, the contributions of 
external current systems have been assumed to be 


zero . 
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THE ARGUS SHELLS 

V/ell defined shells of electrons injected by 
the three high altitude nuclear detonations Argus I, v 
II and III were detected by Satellite 1958e (Explorer IV ) 
/Van Allen, et al., 19527* Each of the observed inter- 
sections of the satellite trajectory with an Argus 
shell should have the same value of L. Figure 4 
shows the L values obtained for the observed inter- 
sections as a function of geographic longitude. The 
average L values for the shells are 1.715* 2.115 and 
2.16 earth radii for Argus I, II and III respectively. 

In addition to the scatter due to measurement in- 
accuracies, systematic variations with longitude of 
approximately +1*5 percent can be discerned. By 
comparison, Pennington /I 961 J has found that the 
trace of the data points in the northern hemisphere 
for Argus III projected upon the surface of the 
earth deviates more than +3*5 degrees from the earth 
trace of the be, 'it possible offset dipole shell. A 
latitude error of + 5*5 degrees for the Argus III 
shell corresponds to an error in L of about + 15 per- 
cent . 

The range of B for the data points on each 
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shell is not large. The deviations of the L values 
in Figure 4 are therefore almost entirely due to 
inaccuracies in the surface field upon which the 
magnetic field representation is based. It is 
interesting to note that the deviations in the L 
values are largest in the regions where the dis- 
crepancies between the various available maps of 
the surface field are also large. No definite time 
dependence can be discerned in the data shown in 
Figure 4. An upper limit on the rate of change in 
L is 0.001 earth radii per day. 
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PROTON INTENSITIES IN THE INNER ZONE 

The most complete set of data pertaining to 
the spatial distribution of particles in the inner 
zone presently available was obtained by the detectors 
on Satellite 1958 Epsilon /y an Allen, Mcllwain and 
Ludwig, 19587* The results of a preliminary analysis 
of the Geiger tube data are shown in Figures 5 and 6. 
The geometric factors and proton thresholds are 

/ ' o 

approximately 0.54 cm and 51 Mev for the unshielded 
counter and 0.62 cm and 45 Mev for the shielded 
counter. Comparison with the scintillation counter 
data indicates that most of the Geiger tube counting 
rates were probably due to penetrating protons rather 
than bremsstrahlung from the very much higher fluxes 
of non-penetrating electrons. 

The data points shown in Figures 5 and 6 were 
obtained from recordings made at stations in North 
and South America, Europe, Africa, Asia, Australia, 

New Zealand and Hawaii. Data obtained during the 
periods immediately following high altitude nuclear 
detonations have been excluded. 

The general physical configuration of the 
lower part of the inner zone is illustrated in 
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Figures 7 and 8 where the B versus L curves in 
Figures 5 and 6 have been transformed into polar 
coordinates in the manner previously described. 

The dashed lines approximately represent the maximum 
excursions of the earth's surface in these coordinates. 

Some of the aspects of these graphs which will 
be subjects for future papers include: 

(a) Determination of angular distributions in 
the manner outlined by E. C. Ray /196 o7* 

(b) The spatial dependence of proton injection 
and loss mechanisms. 

(c) The spatial dependence of the proton 
energy spectrum. 

(d) The time dependence of the Intensities. 

(e) Comparison with data obtained by instru- 
ments on other space vehicles. 
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SOLAR COSMIC RAYS 

The application of the B and L coordinate 
system is by no means limited to trapped particles. 
In general, it can be usefully employed in the study 
of any phenomenon in which motion tends to be along 
lines of force. Some obvious examples are auroral 
particles, whistler mode wave propagation, and low 
energy cosmic rays. 

The latitude dependence of solar cosmic rays 
is of particular interest in that it can be used to 
obtain information about both the energy spectrum of 
the particles and the distortion of the earth 1 s 
magnetic field by external current systems. An 
example of the latitude dependence of solar cosmic 
rays is shown in Figure 9- The left side of this 
figure shows the intensity versus time measured by 
the shielded Geiger counter aboard Satellite 1959 
Iota (Explorer VII ) during a solar cosmic ray event 
and the right side shows the same data points 
plotted versus L. 
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TIME VARYING MAGNETIC FIELD 

The magnetic shells labelled by L values 
greater than 3 earth radii are probably distorted 
by an important amount during magnetic storms. It 
would be difficult if not impossible to calculate 
true instantaneous values of I during a magnetic 
storm. It is suggested that the values calculated 
for magnetically quiet periods be used as a fixed 
reference system in which the perturbations on the 
various phenomena can be studied. 

COMPUTATION OF I 

1/2 

Direct evaluation of I = <J> ( 1 - B^/B) ds 
entails a considerable amount of computation. A 
number of relatively fast methods of computing I 
have been devised. At least one of these methods 
(which in general involve various means of interpo- 
lation) will probably be made available in the form 
of a computer program written in the Fortran language. 
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CONCLUSION 

The choice of B and L as the primary coordinates 
is of course quite arbitrary. One possible alternate 
for B is B/B 0 = L^B/M where B Q = M/L^ is the equatorial 
value of B on the line of force. It should be noted 
that any pair of the quantities B, L, R, Tv and B/B q 
can be used in place of B and L. To avoid confusion 
and to facilitate the comparison of different sets of 
data, it is suggested that a pair of coordinates other 
than B and L be employed only if a substantially 
improved presentation is obtained. 
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Figure 1 

Figure 2 
Figure 3 

Figure 4 
Figure 5 

Figure 6 
Figure 7 


FIGURE CAPTIONS 


The mapping of the polar coordinates R and 
X on to the B,L plane according to the 
transformation 

B = (4 • f ) 1/2 , R = L COB 2 

The variation in L along six lines of 
force. Several of the worst cases are 
included. 

The maximum deviation in L along lines 
forming magnetic shells labelled by 
L = 1.2, 1.5, 2.2 and 3«0 earth radii 
plotted versus the geographic longitude at 
the northern ends of the lines. A positive 
deviation means that L increased toward 
higher values of B. 

The L values for the observed intersections 
of the Explorer IV trajectory with the 
three Argus shells are plotted versus 
geographic longitude. 

Contours of constant true counting rate of 
the unshielded Geiger counter in Explorer IV. 
The points shown correspond to data obtained 
over a wide range of geographic latitude and 
longitude. 

Contours of constant true counting rate of 
the shielded Geiger counter in Explorer IV. 

The contours of constant counting rate 
shown in Figure 5 transformed to the polar 
coordinates R and X. 
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figure 8 


Figure 9 


FIGURE CAPTIONS 
(continued) 


The contours of constant counting rate 
shown in Figure 6 transformed to the polar 
coordinates R and A. 

The left side of this figure shows the 
omnidirectional intensity measured by the 
shielded Geiger counter aboard Explorer VII 
as a function of time during a solar cosmic 
ray event. The right side of the figure 
shows the same data points plotted versus L 
On some occasions the modification of the 
earth's field by external current systems 
permits solar particles to arrive along 
lines of force labelled by considerably 
lower L values. 
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TABLE I 


I ? B/M 

l ? B/M 

T ? B/M 

L ? B/M 

2.37529 x 10" 15 

1.00002 

1.93121 X 10" 5 

1.17553 

1.72697 x io‘ 12 

1.00016 

2.65795 x 10" 5 

1.19618 

3.70298 x 10 -11 

1.00045 

3.60595 x 10 

1.21827 

2.37086 X 10 ” 9 

1.00180 

4.82943 x io"- 5 

1.24187 

2.70390 X 10 " 8 

1.00406 

6.39281 X 10' 5 

1.26706 

1.52215 X 10" 7 

1.00722 

8.37277 x 10 

1.29392 

5.82114 x 10~ 7 

1.01131 

1.08602 x 10 " 2 

1.32252 

1.74352 x 10" 6 

1.01632 

1.39624 x 10 * 2 

1.35298 

4.41257 x 10' 6 

1.02228 

1.78056 X 10 " 2 

1 . 38538 

9.87364 X 10 

1.02919 

2.25380 X 10 ' 2 

1.41985 

2.01130 x 10 

1.03707 

2.83332 X 10 “ 2 

1.45650 

3.80506 X 10" 5 

1.04596 

3.53938 x 10" 2 

1.49547 

6.78113 x 10" 5 

1.05585 

4.39567 x 10 “ 2 

1.53691 

1.15049 x 10 ~ 4 

1.06680 

5.42973 x 10 " 2 

1.58097 

1.87314 X 10 ’ 4 

1.07881 

6.67358 x 10 " 2 

1.62782 

2.94482 x 10 “ 4 

1.09194 

8.16455 x 10 _2 

1.67767 

4.49250 X 10 ” 4 

1.10620 

9.94598 x 10 ' 2 

1.73070 

-4 

6.67702 x 10 

1.12165 

1.20681 x lo’ 1 

1.78716 

9.69961 x 10 " 4 

1.13832 

1.45893 x 10 -1 

1.84729 

1.38095 X 10 " 5 

1.15626 

1.75772 X lo” 1 

1.91135 
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TABLE I 
(continued) 


Jr B/M 


1? B/M 


B/M 


It B/M 


2.11103 X 10" 1 

1.97966 

2.52798 X 10 ” 1 

2.05255 

3.01912 x 10’ 1 

2.13036 

3.59672 x 10" 1 

2.21352 

4.27504 X 10 _1 

2.30245 

5.07058 x 10' 1 

2.39766 

6.00263 x 10 ^ 

2.49967 

7.09352 x lo" 1 

2.60911 

8.36954 x 10 _1 

2.72664 

9.86054 x io" 1 

2.85300 

1.16025 

2.98905 

1.36367 

3.15571 

1.60116 

3.29403 

1.87843 

3.46520 

2.20216 

3*65056 

2.58020 

3.85161 

3.02183 

4.07005 

3.55799 

4.30783 

4.14164 

4.56716 


4.84816 

4.85057 

5.67587 

5.16095 

6.64657 

5.50164 

7*78642, 

5.87648 

9.12680 

6.28994 

1.07055 x 10 

6.74720 

1.25680 x 10 

7.25432 

1.47699 x 10 

7.81840 

1.73787 x 10 

8.44/81 

2.04766 X 10 

9.15245 

2.41655 x 10 

9.94413 

2.85705 x 10 

1.08369 x 10 

3.38471 x 10 

1.18478 x 10 

4.01899 x 10 

1.29975 x 10 

4.78427 x 10 

1.43104 x 10 

5.71146 x 10 

1.58179 x 10 

6.83985 x 10 

1.75576 X 10 

8.21991 x 10 

1.95767 x 10 

9.91695 x 10 

2.19347 x 10 
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TABLE I 
(continued) 


I? B/M 

1.20162 X 10 2 
1.46501 X 10 2 
1.79085 x 10 2 
2.20532 X 10 2 
2.45179 X 10 2 
2.73402 x 10 2 
5.05069 X 10 2 
5.41513 x 10 2 
5.82642 x 10 2 
4.50241 x 10 2 
4.84511 x 10 2 
5.47285 x 10 2 
6.19345 x 10 2 
7.05884 x 10 2 
a. 01427 x 10 2 
9.16824 x 10 2 
1.05125 x 10^ 

1.21182 X 10 5 
1.40095 x io- 5 


l? b/m 

2.47069 x 10 
2.79896 x 10 
5.19079 x 10 

5.66255 x 10 
3.93495 x 10 
4.25601 x 10 
4.56967 x 10 
4.94049 x 10 
5.35387 x 10 
5.81615 x 10 

6.53486 x 10 
6.91896 x 10 
7.57921 X 10 
8.52855 x 10 
9.18267 x 10 

1.01607 x 10 2 
1.12862 x 10 2 
1.25885 x 10 2 
1.41052 x 10 2 


I 5 B/M 

1.62942 X 10 5 

1.90201 X 10 5 
2.25574 X 10 5 
2.64005 X. 10* 
5.14500 X 10 5 

5.76561. X 10 5 
4.55081, X 10 5 
5.54426 X 10 5 , 
6.85524 x 10 5 
8.51087 x 10 5 
1.07567 x 10 4 
1.37789 x 10 4 

1.79974 x 10 4 
2.39555 X 10 4 
5.27489 X 10 4 
4.60402 X 10 4 

6.73169 x 10 4 

1.02942 X 10 5 

1.67841 x 10 5 


1? B/M 

1.58767 x 10 2 
1.79669 x 10 2 
2.04485 x 10 2 
2.54175 x 10 2 
2.70014 X 10 2 

5.15689 x 10 2 
5.67477 x 10 2 

4.54504 X 10 2 

5.19133 x 10 2 
6.27578 x 10 2 
7.68892 x 10 2 
9.56605 x 10 2 
1.21155 x 10 5 
1.56685 x 10 5 
2.07745 x 10 5 
2.85877 x 10 5 
4.02645 x 10 5 
5.98740 x 10 5 
9.47159 x 10 5 
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TABLE I 
(continued) 


I 5 B/M 

1 ? B/M 

I 5 B/M 

1 ? B/M 

2.96984 X 10 5 

1.65048 x 10 ^ 

4.09266 X 10 7 

2.00752 x 10 6 

5.95512 X 10 5 

5.17247 x 10 4 

1.91226 x 10 8 

9.28015 x 10 6 

1.44455 x 10 6 

7.49155 x 10 ^ 

6.50554 x 10 8 

5.12969 x 10 7 

5.00559 X 10 6 
1.47748 x 10 7 

2.51887 X 10 5 
7.52706 X 10 5 

5.20779 X 10 9 

2.50188 x 10 8 
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TABLE II 



X < -16 

-16< X<0 

0< X<8 

8< X421 

X >21 

a o 

0.294 

0.62290 

0.62291 

1.0824 

-5.04 

a l 

0.550 

0.45551 

0.45416 

0.20595 

1.00 

a 2 

0 

1.4495 X lo' 2 

1.5680 X 10 " 2 

5.4145 x 10 " 2 

0 

a 5 

0 

1.2154 X 10 

1.4784 x 10 

-9.5218 x 10" 4 

0 

a 4 

0 

5.9474 X 10 " 5 

1.2415 x 10' 5 

-5.6851 x 10 " 5 

0 

a 5 

0 

i.5567 X 10 " 6 

-8.1278 x 10 “ 6 

2.7879 x 10 ' 6 

0 

a 6 

0 

1.5845 x 10 “ 8 

1.46o4 x 10 " 7 

-5.4751 x 10 ” 8 

0 
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